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This  report  describes  an  algorithm  that  is  used  to  determine  the  maximum  likelihood  position  of  an  aircraft 
from  the  measurement  of  one  radar's  range,  azimuth,  and  elevation  and  from  an  altitude  measurement  of  the 
aircraft.  The  variances  of  these  measurements  must  also  be  known. 

The  report  also  shows  how  much  the  altitude-aided  algorithm  improves  the  accuracy  as  compared  to  a  one- 
radar  determination. 

It  was  concluded  that: 

1.  The  position  error  in  the  low  radar  elevation  angle  ACMTS  (Air  Combat  Maneuvering  Test  System)  geometry 
at  the  Pacific  Missile  Range,  obtained  by  using  a  telemetered  altitude  measurement  and  the  R,  A,  and  E  from 
one  radar  in  the  algorithm  described  in  the  report,  could  be  reduced  from  1/2  to  1/S  that  of  the  position 
error  obtained  when  only  the  R,  A,  E  measurement  from  one  radar  is  used. 

2.  The  altitude  measurement  greatly  reduces  the  error  in  determination  of  the  position  point  when  the  aircraft 
being  tracked  is  a  large  distance  from  the  radar. 


DD  .”*“..1473  <PA6E*> 


S/N  0t01.fl07.680l 


UNCLASSIFIED 

Security  Classification 


UNCLASSIFIED 

Stcutlty  CUnUlctliorr 


KEY  WORDS 


Altitude  measurement 

Low  elevation  angle  radar  tracking  improvement 

'Maximum  likelihood  estimate 

Multivariate  normal  probability  density  function 

Radar  measurement 

Random  vector 

Real  time  application 

Tracking 

Variance -co variance  matrix 


DD  ,’r..1473  <«°» 

(PAGE  £) 


. .  _ 


CONTENTS 

Page 

SUMMARY . 1 

INTRODUCTION . 3 

METHOD . 4 

EFFECTS  OF  THE  ALGORITHM  AND  IMPROVEMENT  OF  TRACKING  DATA  ....  9 

CONCLUSIONS . 10 

REFERENCES . ;  10 

APPENDIX 

Flowcharts . 25 


|.  FIGURES 

1 .  Example  of  a  Geometric  Situation  Used  to  Determine  Algorithm 

t.  Effects  and  Improvement  in  Accuracy . 11 

f  2.  Plots  Showing  the  Effects  of  the  Algorithm  With  Respect  to  Variances  , 

in  the  Elevation  Measurement . 12 

l  3.  Plot  Showing  the  Increase  in  Precision  of  the  Algorithm  Over  the 

}  Radar-Only  Solution . 15 

‘  4.  Plot  Showing  the  Increases  in  Position  Accuracy  as  Range  Increases . 16 

5.  Plots  Showing  the  Cartesian  Coordinate  That  the  Altitude  Measure¬ 
ment  Aids  .  17 

6.  Example  of  a  Geometric  Situation  Used  to  Determine  Improvement 

in  Algorithm  Accuracy . 23 


PACIFIC  MISSILE  RANGE 

Point  Mugu,  California 


ALTITUDE-AIDED  RADAR  TRACKING 

By 

L.  L.  GOERTZEN 


SUMMARY 

This  report  describes  an  algorithm  that  is  used  to  determine  the  maximum  likelihood  posi¬ 
tion  of  an  aircraft  from  the  mea:urement  of  one  radar’s  range,  azimuth,  and  elevation  and  from 
an  altitude  measurement  of  the  aircraft.  The  variances  of  these  measurements  must  also  be 
known. 

The  report  also  shows  how  much  the  altitude-aided  algorithm  improves  the  accuracy  as 
compared  to  a  one-rada”  determination. 

It  was  concluded  that: 

1 .  The  position  error  in  the  low  radar  elevation  angle  ACMTS  (Air  Combat  Maneuvering  Test 
System)  geometry  at  the  Pacific  Missile  Range,  obtained  by  using  a  telemetered  altitude 
measurement  and  the  R,  A,  and  E  from  one  radar  in  the  algorithm  described  in  the  report, 
could  be  leduced  from  1/2  to  1/5  that  of  the  position  error  obtained  when  only  the  R,  A, 
E  measurement  from  one  radar  is  used. 

2.  The  altitude  measurement  greatly  reduces  the  error  in  determination  of  the  position  point 
when  the  aircraft  being  tracked  is  a  large  distance  from  the  radar. 
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INTRODUCTION 

This  report  describes  an  algorithm  that  is  ,«s2d  to  determine  the  maximum  likelihood  posi¬ 
tion  of  an  aircraft  from  the  measurements  of  ran^e,  azimuth,  and  elevation  from  one  radar  and 
from  the  altitude  measurement  of  the  aircraft.  This  algorithm,  a  direct  application  of  informa¬ 
tion  in  reference  1 ,  would  be  easy  to  extend  to  include  data  from  several  radars.  However,  this 
report  is  written  with  the  ACMTS  (Air  Combat  Maneuvering  Test  System)  in  mind.  That  system 
requires  that  from  two  to  six  aircraft  be  tracked  simultaneously;  hence  it  appears  that  the  situa¬ 
tion  where  more  than  one  radar  will  track  a  single  aircraft  will  seldom  be  realized. 

In  reference  2,  recommendation  is  made  to  use  a  radar  altitude-aided  position  determina¬ 
tion  algorithm  that  neglects  the  elevation  measurement  completely.  Reference  2  rejected  the 
method  of  weighting  all  the  data  according  to  anticipated  accuracies  because  that  type  of  method 
complicated  the  calculations  and  did  not  appear  to  offer  a  significant  improvement  in  the  accu¬ 
racy  of  low-altitude  ,neasuu.ments  over  the  method  which  neglected  the  elevation  measurement. 

The  maximum  likelihood  estimate  may  require  more  calculations  to  accomplish;  specifi¬ 
cally  several  more  sine,  cosme  functions  must  be  calculated.  However,  this  procedure  is  not 
complicated  and  certainly  is  usable  in  real-time.  The  procedure  outlined  in  this  report  has  the 
following  advantages  over  the  technique  which  neglects  the  elevation  angle  measurement: 

1.  The  elevation  angle  (even  at  small  elevations)  may  be  a  better  measurement  than  the  alti¬ 
tude  measurement.  This  depends  upon  the  location  of  the  aircraft  with  respect  to  the 
radar.  For  example,  when  the  aircraft  is  less  than  15  miles  from  the  radar,  the  elevation 
angle  measurement  v,  >th  a  standard  deviation  of  0.5  degree  is  better  than  an  altitude  meas¬ 
urement  with  a  standard  deviation  of  100  feet. 

2.  It  is  a  good  general  purpose  procedure.  It  can  use  combined  altitude  and  elevation  data 
when  these  measurements  are  of  nearly  equal  accuracy. 

3.  It  allows  a  continuous  transition  from  inaccurate  elevation  angle  data  to  accurate  elevation 
angle  data.  When  the  accuracy  of  the  elevation  angle  gains  quality  and  should  be  used  in 
the  position  solution,  this  algorithm  does  not  have  to  be  disregarded  and  a  new  procedure 
adopted.  This  algorithm  provides  a  way  to  use  accuracy  evaluations  of  the  elevation  3ngle 
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data  as  an  operation  proceeds  and  to  weight  the  elevation  data  accordingly  in  Lie  position 
solution. 

Reference  2  and  the  paragraphs  of  this  report  that  describe  the  effects  of  the  algorithm 
indicate  that  altitude-aided  radar  tracking  greatly  increases  the  accuracy  over  tracking  that  uses 
radar  data  alone.  With  the  ACMTS,  where  altitude  data  from  a  telemetry  source  is  available  in 
the  computer,  it  would  be  negligence  if  an  altitude-aided  radar  tracking  procedure  were  not 
included  in  the  ACMTS  software. 

This  report  will  proceed  by  first  mathematically  describing  the  algorithm;  then  the  effects 
and  improvements  resulting  from  the  algorithm  will  be  shown  by  examining  the  following: 

1 .  The  influence  of  the  altitude  accuracy  and  altitude  measurement  in  the  maximum  likeli¬ 
hood  range,  azimuth,  and  “levation  as  compared  to  the  range,  azimuth,  and  elevation  of 
the  radar  alone. 

2.  The  influence  of  the  altitude  accuracy  on  the  accuracy  of  the  maximum  likelihood  posi¬ 
tion  estimate. 

3.  The  determination  of  the  slant  ranges  at  which  altitude-aided  tracking  gives  significant 
increases  in  the  accuracy  of  the  position  point. 


METHOD 


Suppose  that  one  radar  measures  the  position  of  an  aircraft  and  that  the  altitude  of  the 
object  above  the  earth  is  measured.  It  is  desired  to  combine  these  measurements  into  an  optimal 
position  determination  of  the  object. 


It  is  assumed  that 

1 .  The  radar  measurement  is  independent  of  the  altitude  measurement,  and 

2.  The  probability  density  functions  of  the  radar  (fR)  and  altitude  (fR)  measurements  are 
gaussian,  i.e., 
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where 


(R,  A,  E)  is  the  .andom  vector  representing  the  radar  measurements 
Mr  is  the  variance-covariance  matrix  of  (R,  A,  E) 

(Mr./wa.  /*e)  is  the  mean  of  (R,  A,  E) 

(x,y,z)  is  the  random  vector  ?epresenting  the  altitude  measurements 
z  is  the  altitude  measurement 

Since  x  and  y  are  not  measured,  their  variances  are  infinite 
Mh  is  the  variance-covariance  matrix  of  (x,y,z) 

(A*x,  /*y,  A*r)  «  the  mean  of  (x.y.z) 

3.  The  spherical  variance-covariance  matrix,  MR,  can  be  transformed  to  a  Cartesian  variance 
matrix  by  the  linear  map  approximation  of  the  first  order  Taylor  Series. 

4.  The  latitude  and  longitude  of  the  location  of  the  altitude  measurement  are  known.  This 

is  a  reasonable  assumption  since  (R,  A,  E)  measurements  can  be  used  to  calculate  this  point 
and  the  variance  in  the  altitude  will  be  insensitive  to  small  errors  in  the  location  of  the 
altitude  measurement. 

mr  can  be  propagated  to  a  Cartesian  variance-covariance  matrix,  MRH,  referenced  to  the 
point  on  the  earth  where  the  altitude  was  measured,  i.e., 
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RH 
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where  T2  and  T3  are  3x3  matrices  (see  reference  3),  and  T*  and  T^  are  the  transpositions  of 
the  matrices  T2and  Tj. 

Tj  (1,  I)  =  cos  E  sin  A 

T3  (1,  2)  -  R  cos  E  cos  A 

Tj  (1,  3)  =  -R  sin  E  sin  A 

Tj  (2,  1)  =  cos  E  cos  A 

T3  (2,  2)  =  ~R  cos  E  sin  A 

Tj  (2,  3)  =  -R  sin  E  cos  A 

Tj  (3,  1)  =  sin  E 


T3  (3,  2)  =  0 

T3  (3,  3)  =  R  cos  li 

T2(l,  1)  =  cos  (Xp  -  Xr) 

T2  (1,  2)  =  sin  0r  sin  (Xp  -  Xf) 

T2  (1,  3)  =  -  cos  <pt  sin  (Xp  -  X,j 

T2  (2,  1)  =  -  sin  0psin(Xp-Xf) 

T2  (2,  2)  =  cos  4>t  cos  0 p  +  sin  0f  sin  0p  cos  ^Xp  -  \  j 

T2  (2,  3)  =  sin  0t  cos  0p  -  cos  0f  sin  0p  cos  |Xp  -  Xfj 

T2  (3,  l)  =  cos0psin(Xp-Xt) 

T2  (3,  2)  =  cos  0r  sin  0p  -  sin  0f  cos  0p  cos  (  Xp  -  Xf  j 

T2  (3,  3)  =  sin  0f  sin  0p  +  cos  0f  cos  0p  cos  (Xp  -  Xfj 

where  0r  and  Xf  are  the  latitude  and  longitude  of  the  origin  of  the  radar  and  0p  and  Xp  are  the 
latitude  and  longitude  of  the  position  of  the  aircraft. 

The  probability  density  function  of  the  radar  measurements  in  the  Cartesian  coordinate 
system  whose  orty'a  is  at  the  point  on  the  earth  where  the  altitude  was  measured  is 


fRC  (XR’yR’ZR) 
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From  assumption  1*  the  probability  density  function  of  x,  y,  z,  xR,  yR,  zR  is 


xR"^x 
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Radar  measurement  is  independent  of  the  altitude  measurement. 
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As  in  reference  1,  pages  4  through  7,  a  mean  0*x,  My,  Pz)  is  found  such  that  for  the  prob¬ 
ability  density  function  f,  the  data  point  (x,  y,  z,  xR,  yR,  zR)  is  most  likely  to  approximate 
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M^1  becomes 
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since  x  and  y  are  not  measured;  their  variances  are  assumed  to  be  infinite. 


Now  directly  from  reference  1 


c6  =  S23 

C7  ~  SH  XR  +  S12  yR  +  S13  ZR 
c8  =  S12  XR  +  S22  yR  +  S23  ZR 


c9  =  S13  XR  +  S23  yR  +  S33  ZR  + 

z 

The  variance-covariance  mat:ix  of  (a,  fi  t  u)  is 

x  y  z 


Fee  the  appendix  for  logical  flow  and  arrangement  of  the  equations  presented  into  a  com¬ 
puter  program. 
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EFFECTS  OF  THE  ALGORITHM  AND  IMPROVEMENT  OF,  TRACKING  DATA 

The  geometrical  situation  diagramed  in  figure  1  was  used  in  the  generation  of  figures  2 
and  3.  It  is  assumed  that  the  radar  measures  an  aircraft  at  (R,  A,  E)  =  (SO  nautical  miles,  20 
i  degrees,  10  degrees)  when  the  altitude  measurement  on  the  aircraft  is  1  nautical  mile.  These 

!  two  measurements  were  made  inconsistent  purposely  so  that  it  can  be  observed  how  the  algo¬ 

rithm  solution  converges  to  the  radar  solution  as  the  variance  in  the  altitude  measurement 
increases.  If  the  radar  were  making  a  consistent  measurement  to  the  altitude  of  1  nautical  mile, 
the  elevation  angle,  E,  would  be  approximately  0.7  degree. 

Figure  2  compares  the  algorithm’s  determination  of  R,  A,  and  E  respectively  to  the  R,  A, 
and  E  of  radar  measurements.  Figure  2  also  shows  how  the  algorithm  converges  to  the  radar 
measurements  as  the  variance  in  the  altitude  measurement  increases,  as  well  as  how  the  altitude 
measurement  and  its  accuracy  affect  the  range,  azimuth,  and  elevation  estimates.  The  symbols 
1,  2,  and  3  in  figure  2  represent  different  standard  deviations  in  the  elevation  angle:  1  assumes 
a  aE  =  10  mils,  2  assumes  a£  =  1  mil,  and  3  assumes  ri£  =  0.1  mil.  The  standard  deviations  of 

range  and  azimuth  of  the  radar  measurements  are  held  constant  at  10  feet  and  0.1  mil. 

Figure  3  shows  how  position  error*  of  the  algorithm  converges  to  the  position  error  of 
!  radar  measurement  as  the  uncertainty  in  the  altitude  increases.  The  increase  in  accuracy  of  the 

algorithm  with  respect  to  the  radar-only  measurement  as  a  function  of  altitude  variance  is  shown 
also.  The  symbols  1,  2,  and  3  in  figure  3  represent  oE’s  of  10  mils,  1  mil,  and  0.1  mil  where 

aR  equals  10  feet  and  aA  equals  0.1  mil  remain  constant.  The  horizontal  lines  at  position  errors 

of  185.1  feet,  85.4  feet,  and  39.6  feet  indicate  the  position  errors  of  the  radar-only  measure¬ 
ment  for  the  respective  standard  deviations  in  elevations  of  10  mils,  1  mil,  and  0.1  mil. 

Figure  4  shows  how  the  altitude-aided  radar  tracking  algorithm  increases  the  position  accu¬ 
racy  as  compared  to  radar-only  tracking.  The  symbol  1  in  figure  4  describes  the  position  error 
of  the  radar  alone  as  a  function  of  slant  range,  and  the  symbol  2  describes  the  position  error  of 
the  altitude-aided  radar  position  determination  as  a  function  of  slant  range.  The  following  typi¬ 
cal  values  for  standard  deviations  in  range,  azimuth,  elevation,  and  altitude  were  used  in  calcu¬ 
lating  figure  4:  0R  =  20  feet,  oA  =  0.2  mil,  0£  =  0.2  mil,  and  0R  =  20  feet. 

Figure  5  shows  how  the  x,  y,  z  rectangular  radar  components  are  affected  by  the  altitude- 
aided  radar  tracking  computations.  An  assumed  trajectory  of  z;  =  30,000  feet,  x5  =  y}  =  i  (10,000) 

feet;  i  =  1 , 2, . . . ,  200  was  converted  to  range,  azimuth,  elevation,  and  altitude  parameters.  (See 
figure  6.)  Noise  with  a  standard  deviation  of  20  feet  was  placed  on  range,  0.2  mil  was  placed  on 
azimuth,  0.2  mil  was  placed  on  elevation,  and  20  feet  was  placed  on  altitude.  These  noisy  param¬ 
eters  were  input  into  the  altitude-aided  tracking  algorithm  and  the  x,  y,  z  outputs  from  the  algo¬ 
rithm  were  differenced  with  the  nominal  x,  y,  z  and  plotted  versus  i  =  1, 2, . . . ,  200.  The  resulting 
plots  are  shown  in  figure  5,  parts  a,  c,  and  e.  To  obtain  comparisons  to  the  radar-only  position 
determination,  the  noisy  range,  azimuth,  and  elevation  parameters  were  used  to  compute  x,  y,  z 


*  Position  error  is  defined  as  the  radius  of  a  sphere  whose  volume  is  that  of  the  volume  of  the  69  percent 
probability  ellipsoid.  For  discussions  o  probability  ellipsoids  see  reference  4. 
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and  these  were  differenced  with  the  nominal  x,  y,  z  and  plotted  versus  i  =  1, 2, . . .  ,200.  The 
resulting  plots  are  shown  in  figure  5,  parts  b,  d,  and  f.  Comparing  the  plots  in  figure  5  shows  that 
the  altitude-aided  tracking  technique  has  little  affect  on  the  precision  of  x  and  y  but  greatly 
increases  the  precision  of  z. 


CONCLUSIONS 

1.  The  position  error  in  the  low  radar  elevation  angle  ACMTS  geometry  at  PMR  obtained  by 
using  a  telemetered  altitude  measurement  and  the  R,  A,  and  E  from  one  radar  in  the  maxi¬ 
mum  likelihood  algorithm  described  in  this  report,  could  be  reduced  from  1/2  to  1/S  that 
of  the  position  error  obtained  when  only  the  R,  A,  and  E  measurement  from  one  radar  is 
used.  The  factors  of  1/2  and  1/S  are  obtained  by  using  figure  3  and  assuming  the  standard 
deviation  in  elevation  is  between  1  to  10  mils  and  the  standard  deviation  in  altitude  is 
between  20  and  100  feet. 

2.  The  altitude  measurement  greatly  reduces  the  error  in  determination  of  the  position  point 
when  the  object  being  tracked  is  a  large  distance  from  the  radar.  (See  figure  4  and  figure 
5  parts  e  and  f.) 

3.  The  position  error  reduction  occurs  in  the  components  in  which  the  altitude  was  measured. 
See  figures  5  and  6. 

4.  Some  altitude-aided  radar  tracking  technique  should  be  incorporated  in  the  ACMTS  real¬ 
time  software.  Two  possibilities  are: 

a.  The  technique  outlined  in  reference  2.  Its  advantage  is  the  shortness  of  computations. 

b.  The  technique  outlined  in  this  report.  Its  advantages  are  that  it  is  more  accurate  than 
the  technique  outlined  in  reference  2  and  is  versatile  in  adjusting  to  changing  accu¬ 
racies  in  the  measured  data. 
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STANDARD  DEVIATION  OF  ALTITUDE  (FEET) 

NOTES:  1.  STANDARD  DEVIATIONS  IN  ELEVATION 
ARE  REPRESENTED  BY  THE  FOLLOWING 
SYMBOLS: 

1  -  -  cte  =  10  MILS 


3  -  -  og  =0.1  MIL 

2.  STANDARD  DEVIATION  IN  RANGE  (oR> 
IS  CONSTANT  AT  10  FEET. 

3.  STANDARD  DEVIATION  IN  AZIMUTH 
(oA>  IS  CONSTANT  AT  0.1  MIL. 


Figure  2.  Plots  Siiowing  the  Effects  of  the  Algorithm  With  Respect  to  Variances  in  the  Elevation  Measurement. 
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STANDARD  DEVIATION  OF  ALTITUDE  (FEET) 


NOTES:  1.  STANDARD  DEVIATIONS  IN  ELEVATION 
ARE  REPRESENTED  BY  THE  FOLLOWING 
SYMBOLS: 

1  --oE  =  :0  MILS 

2  -  -  <’c  =  1  MIL 

3  -  -  og  =0.1  MIL 

2.  STANDARD  DEVIATION  IN  RANGE  (oR) 

IS  CONSTANT  AT  10  FEET. 

3.  STANDARD  DEVIATION  IN  AZIMUTH  (oA) 
IS  CONSTANT  AT  0.1  MIL. 


Figure  2.  Continued. 
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STANDARD  DEVIATION  OF  ALTITUDE  (FEET) 


NOTES: 


1.  STANDARD  DEVIATIONS  IN  ELEVATION 
ARE  REPRESENTED  BY  THE  FOLLOWING 
SYMBOLS: 

1  —  Og  =  10  MILS 

2- -aE  =  1  MIL 

3 - -oF  “0.1  MIL 

2.  STANDARD  DEVIATION  IN  RANGE  (oR) 

IS  CONSTANT  AT  10  FEET. 

3.  STANDARD  DEVIATION  IN  AZIMUTH 
(oA )  IS  CONSTANT  AT  0.1  MIL. 


Figure  2.  Concluded. 
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NOTES:  1.  STANDARD  DEVIATIONS  IN  ELEVATION 
ARE  REPRESENTED  BY  THE  FOLLOWING 
SYMBOLS: 

1  --oE  »  10  MILS 
2--o£  »  1  MIL 
3  -  - oE  =0.1  MIL 

2.  STANDARD  DEVIATION  IN  RANGE  (oR) 

IS  CONSTANT  AT  10  FEET. 

3.  STANDARD  DEVIATION  IN  AZIMUTH 
(oA)  IS  CONSTANT  AT  0.1  MIL. 


Figure  3.  Plot  Showing  the  Increase  in  Precision  of  the  Algorithm  Over  the  Reder-Only  Solution, 
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POSITION  ERROR  (FEET) 


NOISE  X  MINUS  TRUE  X  (FEET! 


NOTES:  1.  MAXIMUM  LIKELIHOOD  ERROR  IS 
REPRESENTED  BY  THE  SYMBOL  M. 

2.  STANDARD  DEVIATIONS  ARE  THE 
FOLLOWING: 

RANGE  - "0  FEET 
AZIMUTH  •  -oA  =0.2  MIL 
ELEVATION  -  -  =0.2  MIL 

ALTITUDE  --aH  =  20  F H ET 

Figure  5.  Plots  Shoeing  the  Cartesian  Coordinate  That  the  Altitude  Measurement  Aids. 


ISE  X  MINUS  TRUE  X  (FEET) 
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X  POSITION  (10,000  FEET) 


NOTES:  1.  RADAR  ERROR  IS  REPRESENTED 
eY  THE  SYMBOL  R. 

2.  STANDARD  DEVIATIONS  ARE  THE 
FOLLOWING: 

RANGE  --oR  =20  FEET 
AZIMUTH --oA  =0.2  MIL 
ELEVATION --oE  =0.2  MIL 


Figure  5.  Continued. 
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NOISE  Y  MINUS  TRUE  Y  (FEET) 


Figure  5.  Continued. 


SE  Y  MINUS  TRUE  Y  (FEET) 


NOISE  Z  MINUS  TRUE  Z  (FEET) 


NOTES:  1.  MAXIMUM  LIKELIHOOD  ERROR  IS 
REPRESENTED  BY  THE  SYMBOL  M. 

2.  STANDARD  DEVIATIONS  ARE  THE 
FOLLOWING: 

RANGE  =20  FEET 
AZIMUTH  --oA  =0.2  MIL 
ELEVATION --<;E  =0.2  MIL 
ALTITUDE --oH  =20  FEET 


Figurt  5.  Continual. 
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NOISE  Z  MINUS  TRUE  Z  (FEET) 


Figure  5.  Concluded. 


Figure  6.  Example  of  a  Geometric  Situation  Used  to  Determine  Improvement  in  Algorithm  Accuracy. 


APPENDIX 

FLOWCHARTS 


Preceding  page  blank 


I  SET  THESE 
I  PARAMETERS 
I  IN  COMMON 
I  BEFORE 
I  ENTERING 
*  THIS 
|  ROUTINE: 


oR  -  STANDARD  DEVIATION  OF  RANGE  ~"| 
aA  .  STANDARD  DEVIATION  OF  AZIMUTH  . 
aE  -  STANDARD  DEVIATION  OF  ELEVATION  • 
aH  -  STANDARD  DEVIATION  OF  ALTITUDE  | 
Ae  -  EARTH  EQUITORIAL  RADIUS  j 

Be  .  EARTH  POLAR  RADIUS  # 

0R  -  RADAR  LATITUDE  1 

4>r  -  RADAR  LONGITUDE  I 

H„  -  RADAR  ALTITUDE  I 


|  AZIMUTH  OF  THE  X-AXIS  OF  ALL  TANGENT  PLANE  I 

I  RECTANGULAR  COORDINATE  SYSTEMS  IS  ASSUMED  TO  BE  1 
90  DEGREES  FROM  TRUE  NORTH.  | 


INPUT: 

R  -  RANGE 
A  .  AZIMUTH 
E  -  ELEVATION 
H  -  ALTITUDE 


CONVERT  FROM  POLAR 
TO  RECTANGULAR 
COORDINATES: 

xR  -  R  COS  E  SIN  A 
yR  -  R  COS  E  COS  A 
zR  -  R  SIN  E 


^R1  ^R'  Hr  AND  X,,,  yR, 
zR  TO  COMPUTE  THE 
GEODETIC  LATITUDE  -  0V 
OF  xR,  yR,  zRt  AND  THE 
LONGITUDE OF  Xjj, 
yR1 rR. 


26 


-  MATRIX  OF  PARTIALS  WHICH  CONVERTS  A  POLAR 
COVARIANCE  MATRIX  TO  A  RECTANGULAR  ONE 

-  COS  E  SIN  A 

■  R  COS  E  COS  A 

-  -R  SIN  E  SIN  A 

-  COS  E  COS  A 

•  -R  COS  E  SIN  A 

-  -R  SIN  E  COS  A 

■  SIN  E 

.0 

-  R  COS  E 


.  MATRIX  WHICH  CONVERTS  THE  RECTANGULAR  COVARIANCE 
MATRIX  IN  THE  RADAR  SYSTEM  TO  A  RECTANGULAR 
COVARIANCE  MATRIX  IN  TANGENT  PLANE  SYSTEM  LOCATED 
AT  ( 0p ,  0p,  0) 

>  COS  (tip  -  tig ) 

-  SIN  0R  SIN  (<ip  -  tig) 

.  -COS  0R  SIN  (tip  -  t iR) 

.  -SIN  0P  SIN  (ip  -  tip) 

.  COS  0R  COS  0P  +  SIN  0R  SIN  0p  COS  (tip  -  tip) 

•  SIN  0R  COS  Op  -  COS  Op  SIN  0p  COS  (0p  -  tiR) 

.  COS  Op  SIN  (tip  -  <t>p) 

.  COS  Op  SIN  Op  -  SIN  Op  COS  0p  COS  (0 P  -  tip) 

.  SIN  Op  SIN  Op  +  COS  Op  COS  0p  COS  (<f>p  -  tiR) 


o 


